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We discuss the properties of the hadron-quark mixed phase in com- 
pact stars using a realistic equation of state of hyperonic matter and 
the MIT bag model. We find that the equation of state of the mixed 
phase is similar to that given by the Maxwell construction, but that the 
mixed phase becomes mechanically unstable if the surface tension of 
£3 ■ the interface between the two pure phases is strong enough. The com- 

position of the mixed phase is very different from that of the Maxwell 
construction; in particular, hyperons are completely suppressed. 
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It is commonly believed that hyperons appear in dense nuclear matter 
at baryon densities above 2-3 times normal density, in spite of some uncer- 
tainties about the nucleon-hyperon (NY) and hyperon-hyperon (YY) inter- 
actions. Many theoretical studies have shown that the hadronic equation 
of state (EOS) becomes very soft once hyperons become components of the 
matter [H [2]. As a major consequence, the maximum mass of neutron stars 
(NS) predicted using the hyperonic EOS may remain below the current ob- 
servational values of about 1.5 solar masses [3]. 

Some authors have suggested that this situation might be remedied by 
considering the yet unknown three-body forces (TBF) among hyperons and 
nucleons pi], while other studies have shown that a quark deconfinement 
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phase transition in hyperonic matter renders the EOS sufficiently stiff again 
to allow NS masses consistent with current data [5]. 

However, the appearance of quark matter (QM) poses the problem of 
an accurate theoretical description of the quark phase, which is so far an 
open question, and furthermore of the details of the phase transition be- 
tween hadronic and quark matter. The purpose of this letter is the study of 
the latter problem, combining a Brueckner-Hartree-Fock (BHF) EOS of hy- 
peronic hadronic matter with the standard phenomenological MIT model for 
the quark phase. In the simplest scenario, the Maxwell construction (MC), 
a sharp transition takes place between the two charge-neutral hadron and 
quark phases, whereas the more general Gibbs (Glendenning) construction 
(GC) [6] allows a mixed phase (MP) containing individually charged hadron 
and quark fractions with various geometrical structures. However, in the 
latter case, electromagnetic and surface contributions to the energy of the 
MP are usually neglected, but could have important effects [3 El [9]. The 
quantitative analysis of these corrections is the purpose of this letter. 

Our theoretical framework for the hadronic matter is the nonrelativistic 
BHF approach [21 [TO] based on microscopic NN, NY, and YY potentials that 
are fitted to scattering phase shifts, where possible. Nucleonic three-body 
forces (TBF) are included in order to (slightly) shift the saturation point 
of purely nucleonic matter to the empirical value. It has been demonstrated 
that the theoretical basis of the BHF method, the hole-line expansion, is well 
founded: the nuclear EOS can be calculated with good accuracy in the BHF 
two hole-line approximation with the continuous choice for the single-particle 
potential, since the results in this scheme are quite close to the full convergent 
calculations which include also the three hole-line contributions [TO], [XT] . Due 
to these facts, combined with the absence of adjustable parameters, the BHF 
model is a reliable and well-controlled theoretical approach for the study of 
dense baryonic matter. 

The basic input quantities in the Bethe-Goldstone equation are the NN, 
NY, and YY potentials. In this work we use the Argonne Vis NN potential 
[T2] supplemented by the Urbana UIX nucleonic TBF of Refs. [T3] and the 
Nijmegen soft-core NSC89 NY potentials [H] that are well adapted to the 
existing experimental NY scattering data and also compatible with A hy- 
pernuclear levels [T5] [16]. With these potentials, the various G- matrices are 
evaluated by solving numerically the Bethe-Goldstone equation. Then the 
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total nonrelativistic hadronic energy density, e#, can be evaluated: 
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mj + k 2 /2rrii, where the various single-particle potentials are 
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with Ti{k) 
given by 

U l (k) = £ 

i'=n,p,A,£- 

and are determined self-consistently from the G-matrices. 

For the quark EOS, we use the MIT bag model with massless u and d 
quarks and massive s quark with m s = 150 MeV. The quark matter energy 
density can be expressed as a sum of the kinetic term and the leading-order 
one-gluon-exchange term [T7J [18] for the interaction energy proportional to 
the QCD fine structure constant a s , 
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where is the f = u,d,s current quark mass, Xf = kp/rrif, the baryon 

density of / quarks is pf — k^P /3ir 2 , and the bag constant B is the energy 
density difference between the perturbative vacuum and the true vacuum. 

This is clearly an oversimplified model of QM, which will be used in this 
letter to study the generic qualitative features of the hadron-quark phase 
transition in NS matter. In future work we will compare in more detail the 
quantitative results obtained using different, more sophisticated, QM models. 

Figure [1] compares the hadronic BHF EOS and the quark matter EOS 
with different values of the parameters B and a s for beta-stable and charge- 
neutral matter. One can see that the quark EOS approaches that of a rel- 
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ativistic free gas (E/A ~ p# ) with increasing density, while the hyperonic 
EOS is always soft. Consequently the quark deconfinement transition cannot 
occur at too high densities. If we demand the quark and the hyperonic EOS 
to cross, a s should be small and B slightly large, which gives a relatively low 
critical density. Thus the appearance of hyperons is effectively suppressed 
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Figure 1: EOS of hadronic matter (black curves) and of quark matter (col- 
ored curves) with B = 60 MeV/fm 3 (lower curves) and B = 100 MeV/fm 3 
(upper curves) for several values of a s . Hyperons appear at the dotted point 
in hadronic matter. 
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Figure 2: Density profiles and Coulomb potential <fi within a 3D (quark 
droplet) Wigner-Seitz cell of the MP at ps = 0.4 fm -3 . The cell radius and 
the droplet radius are Rw = 26.7 fm and R = 17.3 fm, respectively. 
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Figure 3: EOS of the MP (thick curves) in comparison with pure hadron and 
quark phases (thin curves). The upper panel shows the energy per baryon 
E/A and the lower panel the energy difference between mixed and hadron 
(Pb < 0.44 fm~ 3 ) or quark (p# > 0.44 fm -3 ) phases. Different segments of 
the MP are chosen by minimizing the energy. 
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due to a quark deconfinement transition. In this letter we choose a s = and 
B = 100 MeV/fm 3 . 

The numerical procedure to determine the EOS and the geometrical struc- 
ture of the MP is similar to that explained in detail in Refs. [9] . We employ a 
Wigner-Seitz approximation in which the whole space is divided into equiv- 
alent Wigner-Seitz cells with a given geometrical symmetry, sphere for three 
dimension (3D), cylinder for 2D, and slab for ID. A lump portion made of 
one phase is embedded in the other phase and thus the quark and hadron 
phases are separated in each cell. A sharp boundary is assumed between 
the two phases and the surface energy is taken into account in terms of a 
surface-tension parameter a. We use the Thomas-Fermi approximation for 
the density profiles of hadrons and quarks, while the Poisson equation for the 
Coulomb potential is explicitly solved. The energy density of the mixed 
phase is thus written as 
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where the volume of the Wigner-Seitz cell Vw is the sum of those of hadron 
and quark phases Vh and Vq, and S the quark-hadron interface area. €l 
indicates the kinetic energy density of lepton (only the electron in this work). 
The energy densities e#, £q and are r-dependent since they are functions 
of local densities p a (r) (a = n,p, A, S~, u, d, s, e). For a given density p B , the 
optimum dimensionality of the cell, the cell size Rw, the lump size R, and 
the density profile of each component are searched for to give the minimum 
energy density. The structure of the MP changes from quark droplet to quark 
slab to hadron tube to hadron bubble with increasing baryon density. 

The surface tension of the hadron-quark interface is poorly known, but 
some theoretical estimates based on the MIT bag model for strangelets [17] 
and lattice gauge simulations at finite temperature [19] suggest a range of 
a w 10-100 MeV/fm 2 . We show results using a = 50 MeV/fm 2 in the present 
letter, and discuss the effects of its variation. 

Figure [2] illustrates the outcome of this procedure, showing the density 
profile in a 3D cell for ps = 0.4 fm -3 . One can see the non- uniform den- 
sity distribution of each particle species together with the finite Coulomb 
potential. The quark phase is negatively charged, so that d and s quarks are 
repelled to the phase boundary, while u quarks gather at the center. The 
protons in the hadron phase are attracted by the negatively charged quark 
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Figure 4: Droplet radius (R) dependence of the energy per baryon for fixed 
baryon density ps = 0.35 fm~ 3 and different surface tensions. The quark 
volume fraction (R/Ry^) 3 is fixed for each curve. Dots on the curves show 
the local energy minima. The black line shows the energy of the MC case. 
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Figure 5: Particle fractions of quark and hadron species in the pure hadron 
and quark phases (upper panel) and in the MP (lower panel). In the MC 
the phase transition occurs between the pure phases with pB = 0.34 fm~ 3 
(hadrons) and 0.50 fm~ 3 (quarks). 
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phase, while the electrons are repelled. 

Figure [3] (upper panel) compares the resulting energy per baryon of the 
hadron-quark MP with that of the pure hadron and quark phases over the 
relevant range of baryon density. The thick black curve indicates the case 
of the Maxwell construction, while the colored line indicates the MP in its 
various geometric realizations, starting at ps = 0.326 fm -3 with a quark 
droplet structure and ending at ps = 0.666 fm~ 3 with a hadron bubble 
structure. The energy of the MP is only slightly lower than that of the MC, 
and the resultant EOS is similar to the MC one. However, the structure and 
the composition of the MP are very different from those of the MC case, as 
discussed later. 

If one uses a smaller surface tension parameter a, the energy gets lower 
and the density range of the MP gets wider. The limit of a = leads to a bulk 
application of the Gibbs conditions without the Coulomb and surface effects, 
i.e., the so-called Glendenning construction [6]. On the other hand, using a 
larger value of a, the geometrical structures increase in size and the EOS gets 
closer to that of the MC case. Above a limiting value of a ~ 65 MeV/fm 2 the 
structure of the MP becomes mechanically unstable [5] : for a fixed volume 
fraction (R/Rw) 3 the optimal values of R and Rw go to infinity and local 
charge neutrality is recovered in the MP, where the energy density equals 
that of the MC (see Fig. SJ). 

This mechanical instability is due to the charge screening effect: The 
optimal values of R and Rw are basically determined by the balance between 
the Coulomb energy (~ R 2 in the 3D case) and the surface energy (~ R^ 1 ). 
However, if the charge screening is taken into account, the contribution of 
the screened Coulomb potential <p is strongly reduced when R, Rw —> oo. 
A careful analysis by Voskresensky et al. showed that the Coulomb energy 
changes its sign and behaves like as R — > oo due to the charge screening 
effect [8]. Thus the surface and the Coulomb energy give a local minimum 
below o m 65 MeV/fm 2 , which disappears when the surface energy becomes 
greater than the Coulomb energy above a ~ 65 MeV/fm 2 . This is in contrast 
to the work of Heiselberg et al. [7], neglecting the charge screening effect, 
where there is always a local energy minimum at finite R. The importance of 
the charge screening effect has been also shown in the stability of strangelets 

One notes in Fig. [2] that no hyperons appear in the MP although the 
mean baryon density Pb = 0.4 fm -3 is higher than the threshold density for 
hyperons in pure nucleon matter pg = 0.34 fm -3 (see the black dot in Fig. [I]). 
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In the upper panel, the case of the pure hadron and quark phases is shown 
for comparison. One can see that the particle fractions are very different in 
both cases, in particular a relevant hyperon fraction is only present in the 
hadronic part of the MC. 

Thus we conclude that due to the relatively small magnitudes of the sur- 
face and Coulomb energies, the EOS of the MP is similar to the MC one, but 
the internal structure of the MP is very different. In particular the role of 
hyperons is strongly reduced when we consider the deconfinement transition 
in hyperonic matter. Above a maximum value of the surface tension param- 
eter, the MC is recovered as the physical one, however. These results should 
be important for physical processes like neutrino propagation and baryonic 
superfluidity, besides the maximum mass problem, which will be studied in 
an extended article. 

This work is partially supported by the Grant-in- Aid for the 21st Century 
COE "Center for the Diversity and Universality in Physics" and the Grant- 
in-Aid for Scientific Research Fund of the Ministry of Education, Culture, 
Sports, Science and Technology of Japan (13640282, 16540246). 
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